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ABSTRACT 


The  sound  scattering  due  to  an  ambient  noise  field,  approx¬ 
imated  by  a  squared  cosine  function,  is  considered  for  infinite 
rigid  and  elastic  cylinders  and  rigid  spheres.  For  the  cylinders,  it 
is  assumed  that  the  acoustic  wave  front  is  parallel  to  the  axis  of 
the  cylinder  (normally  incident).  For  this  2issumption,  a  closed 
form  expression  for  the  scattered  sound  field-to-incident  ambient 
noise  field  ( signal- to-noise)  ratio  is  obtained  not  only  for  the  co¬ 
sine  squared  directivity,  but  for  any  arbitraxj'^  directivity  which 
can  be  expressed  in  terms  of  a  Fourier  series.  For  the  sphere,  it 
is  assumed  that  the  noise  is  circumferentially  symmetric  which 
leads  to  a  closed  form  expression  for  the  signal- to- noise  ratio  due 
to  a  cosine  squared  directivity. 


ADMINISTRATIVE  INFORMATION 

This  work  was  carried  out  under  joint  fimding  from  the  Acoustic  Measurement  Facility 
Improvement  Program  (AMFIP)  and  the  Office  of  Naval  Research  (ONR)  Exploratory  Devel¬ 
opment  Program.  The  work  was  carried  out  at  the  Carderock  Division  of  the  Naval  Surface 
Warfare  Center  during  January  1992  to  September  1992. 

INTRODUCTION 

The  ambient  noise  of  the  sea^  is  generated  by  steady  noise  sources,  such  as  surface  winds, 
wave  interactions,  and  distant  ships  as  they  transit  shipping  lanes,  and  by  transient  sources, 
such  as  rain  and  the  calls  of  mzirine  animals.  Although  steady  sources  persist  for  extended 
periods  of  time,  their  variability  leads  to  the  randomness  of  the  ambient  noise.  Intermittent 
or  transient  noise  sources  also  contribute  to  this  variability,  but  cause  greater  uncertainty  in 
expected  noise  levels  since  their  occurrences  axe  less  predictable.  Since  the  ambient  noise  is 
caused  by  a  variety  of  sources,  its  characteristics  change  throughout  the  frequency  spectrum, 
with  different  sources  becoming  dominant  contributors.  Different  characteristics  are  also  ob¬ 
served  at  different  locations.  For  instance,  tidal  currents  are  more  pronounced  in  coastal  waters 
than  in  the  center  of  the  ocean.  In  addition  to  frequency  and  location,  the  noise  characteris- 


1 


tics  are  affected  by  changes  in  sound  transmission  conditions  such  as  those  caused  by  seasonal 
changes. 

Even  though  there  are  a  wide  variety  of  sources  which  can  create  ambient  noise,  the  two 
primary  contributors  are  surface  noises  created  by  winds  and  distant  shipping  noise.  We  con¬ 
sider  frequency  ranges  high  enough  that  contributions  due  to  shipping  noises  can  be  ignored 
since  they  exist  at  much  lower  frequencies.  It  has  been  shown  that  surface  noise  contributions 
extend  over  a  wide  range  of  frequencies  and  on  average  are  circumferentially  symmetric  and 
well  approximated  by  a  cos^  6  distribution.^  Since  the  contributions  from  the  surface  exceed 
those  from  the  sea  bottom,  the  noise  field  will  be  directional. 

The  ambient  noise  field  may  be  thought  of  as  a  source  which,  when  incident  upon  targets 
such  as  cylinders,  will  cause  a  scattered  pressure  field.  There  are  two  primary  concerns  about 
this  scattered  ambient  field.  If  the  scattered  pressure  is  distinguishable  from  the  ambient  noise 
source,  it  may  be  possible  to  use  the  noise  field  as  a  means  of  imaging  the  targets.^  The  other 
concern  is  the  effect  that  the  ambient  noise  scattering  may  have  upon  the  measurement  of  these 
canonical  scatterers  in  scientific  experiments. 

MATHEMATICAL  FORMULATION 

The  mathematical  formulation  of  the  problem  begins  with  a  few  baisic  assumptions;  mea¬ 
surements  are  made  in  a  small  frequency  remge  and  the  linear  principle  of  superposition  can 
be  applied.  Keeping  these  criteria  in  mind,  the  pressure,  ,  at  position  R  due  to  an  ambient 
noise  distribution  can  be  expressed  as  a  superposition  of  plane  waves,  P,  exp[ik-R].  Each  plane 
wave  has  its  own  propagation  direction  and  amplitude  Pi,  but  they  share  the  same  wavenumber 
k.  Therefore,  P;v(R)  can  be  expressed  Jis 

p^(R)  =  j  P^(n')exp[ik'-R]dfi',  (1) 

where  is  the  pressure  amplitude  with  angtilar  dependence  expressed  in  the  appropriate 
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coordinate  system,  k'  is  the  vector  wavenumber  in  the  direction  of  propagation,  and  dQ’  repre¬ 
sents  a  solid  angle  integration  for  the  coordinate  system. 

The  scattered  pressure  from  a  surface  excited  by  an  incident  plane  wave  of  amplitude  P, 
may  be  formally  expressed  as  P,G(R;  Q')  where  G  is  the  impulse  response  function  or  Green's 
function  for  the  siurface.  Applying  the  superposition  principle,  the  scattered  pressure,  ,  at  R 
due  to  an  ambient  noise  field  is  cast  in  the  form 

p,(R)=  (2 


The  effects  of  the  ambient  noise  can  be  determined  by  measuring  the  signal- to-noise  ratio. 
5/iV,  defined  by 


S/N  = 


P(R)-p,(R) 

Pn{^) 


(3) 


where 


P(R)  =  p„(R)-l-p,(R)+p<,(R). 


(4) 


is  the  pressure  detected  by  the  hydrophone  at  R  in  the  presence  of  the  target,  is  the 

pressure  detected  by  the  hydrophone  in  the  absence  of  the  target,  PsCR)  is  the  scattered  am¬ 
bient  pressure  due  to  the  target,  and  Po(R)  is  any  other  observable  pressme.  The  measured 
signal-to-noise  ratio  is  then 


Psm±Po(B  rv 

p.(R)  ■  ^  ^ 

If  the  contribution  of  the  scattered  noise  is  large  enough,  the  ambient  field  can  be  used  for  pas¬ 
sive  imaging.  Far  from  the  scattering  surface  (i.e.  kR  »  1),  the  scattered  pressure  will  dimin¬ 
ish  significantly  as  a  result  of  geometrical  spreading,  and  the  other  observable  pressures  can  be 
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well  estimated  by  using  Eq.  5.  Yet  if  measurements  are  made  close  to  the  scatterer.  the  scat¬ 
tered  noise  and  other  observable  noise  may  not  be  distinguishable  so  that  the  measurements 
may  not  be  statistically  significant. 

Since  the  signaJ-to-noise  ratio  is  a  quantity  of  interest,  we  square  the  expressions  for  p,v(R) 
and  ^^(R)  to  obtain 


ip»(ii)P  =  JJ 


P^{n')P*{n")exp{ik'-K]exp[-ik"-R]dn'dn", 


(6) 


and 


,P.(R)r  = 

2  2 

Due  to  the  randomness  of  the  noise  field,  the  mean  square  pressures,  {|pjv(R)l  )  a^d  (|p5(R)|  ). 
hold  more  physical  meaning  than  a  single  measurement  of  (R)-  Applying  this 

averaging  results  in 

<b„(R)P)=  jJ{PA^')P:(n"))exp[ik'-R]exp[-ik".K]da'dQ'\  (8) 

and 

(|p»(R)l“>  =  j^(P.,(fi')p;(fi"))G(R:n')G-(R;n")<inw.  (9) 

The  term  {P^{Q')P*{Q"))  can  be  expressed  as  a  directivity  function  P(n',f2")  by 

(P^  {n')p;(fi"))  =  lP,f  !>(!)',  n"),  ( lo) 
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where  Pj  is  a  normalization  factor.  Applying  this  to  Eqs.  8  and  9  results  in  the  signd- to- noise 
ratio 


,ffp(n',n")G(R;n')C-(R;n'-Mflw 

'  '  ‘  //P(f!',f!'')exp(ik'R|eicp[-ik"-R|<if!'<ifi"'  ' 

There  are  two  simple  directivity  functions  which  do  not  represent  true  ambient  noise  fields, 
yet  are  still  instructive.  The  first  of  these  is 

=  6{Q'  -  n,)6in"  -  Q,),  (12) 

which  is  the  directivity  for  one  incident  plane  wave  propagating  in  the  direction.  The  sec¬ 
ond  is 


(13) 


which  indicates  no  preferred  direction  and  represents  an  isotropic  noise  field.  A  more  practical 
assumption  used  throughout  this  paper  is  that  the  noise  field  contains  no  correlation  between 
angles.  This  is  equivalent  to  saying 


(14) 


which  simplifies  the  signal-to-noise  ratio  to 


(S/N)^  = 


fI>(Q')G(R;  n')G\R;  Q’)dQ' 
fV(Q'}dQ> 


(15) 
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CYLINDRICAL  SCATTERERS 


INFINITE  RIGID  CYLINDER 

In  cylindrical  coordinates,  a  wave  traveling  perpendicularly  to  the  z  a:>cis  in  the  direction  o' 
with  araphtude  Pi  can  be  represented  by^’^ 

Pinc{R-,<^  ■,<!>')  =  Piexp[ikRcos{<j)’  —  (p)] 

( 16' 

=  Pi 


where  {R,<p)  is  the  location  of  the  receiver,  k  is  the  wavenumber,  and  Jn  is  the  order  cylin¬ 
drical  Bessel  function.  When  an  infinite  rigid  cylinder  of  radius  a  is  excited  by  a  normally  inci¬ 
dent  plane  wave,  the  scattered  pressure  is  given  by^ 


OO 

p,{R,P,4>')  =  -Pi  Y 

n=— OO 


Hi,{ka) 


Hn{kR)e 


(17) 


where  Hn  is  the  cylindrical  Hankel  function  of  n**  order  and  the  primes  on  the  Bessel  and  Han- 
kel  functions  represent  differentiation  with  respect  to  the  arguments  of  these  functions. 

When  more  than  one  normally  incident  plane  wave  excites  this  cylinder,  the  principle  of 
superposition  czin  be  exploited.  With  this  in  mind,  the  normally  incident  field  representing  a 
spatially  continuous  noise  field  of  amplitude  Pj^{<i>')  is  represented  by 


Pn{4>)  =  J Pa{^') exp[tfcRcos(<^'  -  <{>)]  d9! . 


(18) 


The  average  incident  pressure  squared,  (|R/v(^)P)5  continuous  noise  source  is  then 


/3jr/2  fSir/i 

/  dr(.PA4>')p:(r)) 

-ir/2  J  — 


(19) 


exp[ifci2cos(^'  —  <p)]  exp[— ?fcRcos(^"  —  p)]. 
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In  the  same  manner,  the  scattered  pressure  due  to  this  continuous  noise  sourre  represented 
by 


jpM')  t 


.  Jnika) 
H'nika) 


dQ', 


(20) 


The  average  scattered  pressure  squared,  (IPsC^)!^),  is  subsequently  represented  by 


/3tr/2  ^3ir/2 

di,'  /  d0"(PM’)p:(4''')) 

-7r/2  J  — 


jr/2 

oo 


ns=*-oo 
oo 


(21 : 


Define  T„  by 


j"  n  Jn(feQ) 


(22) 


and  define  the  directivity  function  associated  with  the  noise  field,  by 


(e.(<‘')e;(0  =  W,*")|e,r 


(23) 


Using  these  definitions,  Eqs.  19  and  21  simplify  to 


<^") exp(tA:i2(cos(.^'  -  <6)  -  cos(^"  -  <6))]  (24) 

|P;(  y-tr/2  y-ir/2 


and 


(|P5(^)I^) 


/3ir/2  ^3jr/2  »  «. 

<(«'/  d4."V{<t.',n  Yi  11  'P.'Pr  e 

-ir/2  y-ir/2  „=r-oor=-oo 


OO  oo 


(25) 


7 


respectively.  The  next  assumption  made  is  that  the  noise  field  is  spatially  uncorrelated  so  that 


<!>")  =  -r). 


(26) 


Applying  this  to  Eq.  24  yields 


(IPn(*^)I  )  „  f  d(f>'T>{(l)')exp[ikRcos{4>'  —  (i>)\exp[—ikRcos{<p'  —  d>)] 

|PJ  J-r/2 

=  /  d<f>'V{4>'). 

J—kI2 


(27) 


Applying  Eq.  26  to  Eq.  25  and  interchanging  the  order  of  the  summations  and  integrations  re¬ 
sults  in 


ie,P 


n=-oor=-oo 


(28) 


INFINITE  ELASTIC  CYLINDRICAL  SHELL 

The  scattered  pressure  from  an  infinite  elastic  cylindrical  shell  excited  by  an  incident  plane 
wave  propagating  in  the  direction  (f>'  is  given  by^ 


Psei^i^) —  ~Pinc  ^ 


.nHn{kR) 


Jnika)  - 


2pc 


{Zn  +  z„)vkaH'„ika) 


,m(^'  -*) 


(29) 


The  modal  structural  impedance  Zn  is  defined  by 


ip,Cp  h  Q*  -  n^l  + 

Da  D2  —  n2 


(30) 
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where  h/a  is  the  thickness  to  radius  ratio,  =  {h/a)^ /12  is  a  nondimensional  thickness  pa¬ 
rameter,  =  uajcp  is  a  nondimensional  frequency  parameter,  Cp  is  the  compressionaJ  wave 
speed,  and  p,  is  the  density  of  the  structure.  The  specific  acoustic  impedance  z„  is  defined  by 


2n 


=  ipc 


Hnjka) 

Hi,{kay 


where  p  is  the  fluid  density. 

£ 

By  defining  r„  to  be 


(31) 


■MkR) 

HUka) 


Kika)  - 


2pc 

{Zn  +  Zn)TrkaH'^{ka) 


(32) 


Equations  27  and  28  of  the  rigid  cylinder  case  may  be  used  for  an  elastic  cylindrical  shell  by 

E  R  RE 

substituting  T„  for  T„  .  In  the  formulations  to  follow,  r„  will  be  used  to  represent  r„  and  r„  , 
resulting  in  generalized  equations. 


GENERAL  DIRECTIONAL  DIRECTIVITY 

We  will  asstime  that  the  noise  is  normally  incident  on  the  cylinder  so  that  the  directivity 
function  in  cylindriczd  coordinates  is  dependent  only  on  the  angle  <l>.  Although  this  is  not  the 
most  realistic  case,  it  is  a  solvable  problem  and  serves  as  a  first  estimate.  A  general  directional 
directivity  fimction  of  this  sort,  D(^'),  can  be  expressed  in  terms  of  its  Fourier  series  expansion 
as 


_7r/2  <  <^'  <  7r/2 
ff/2  <  ^'  <  37r/2, 


(33) 
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where  Cm  and  bm  are  the  Fourier  series  coefficients  which  are  determined  for  the  ambient  field 
of  interest.  For  this  directivity  function,  the  average  noise  strength  becomes 


/I  |2v  oo  .„/2  oo  -35r/2 

I-PjI  m=-oo  J--K/2  m=-oo  >^w/2 


juid  the  average  scattered  field  yields 


(b.l') 

\P:t 


oo  oo 


=  E  E 


nss  — csQ  rss  — oo 


L  m^  —  oo 


W  .Tfy 

E  -/ 

=  -oo 

oo  - 

^  6„»/ 

m=-oo 


jR  £ 

In  Eq.  35,  T„  represents  either  T„  or  r„  ,  depending  upon  the  case  of  interest. 
Upon  the  simplifications  of  Appendix  B,  these  expressions  become 


l•wl  m=-oo 


odd 


and 


(ip,r) 

1^,1* 


=  )r  E  (<■".+ E 


mss— Qo 

oo 


n=^oo 
oo  oo 


+2  E  E  E 

mss— OO  nss— oorss— oo 

g-«(n-r)^^_jj(n-r+m-l)/2 

(n  —  r  +  m) 


bn  —  r+m,odd  • 


(34) 


(35) 


(36) 


(37) 


Isotropic  Noise  Field 

An  isotropic  noise  field  is  one  in  which  the  noise  is  evenly  distributed  in  all  directions.  The 
directivity,  Pi#o(^'),  is  therefore  a  constant  so  that  the  F^'urier  series  coeflScients  are  oo  =  bo 
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and  Om  =  6m  =  0  Vm  >  1.  The  resulting  isotropic  noise  field  is  given  by 


=  27rao 


(38) 


zuid  the  isotropic  scattered  field  is  given  by 


{\Ps\\so 

\PA" 


oo 

=  27rao 

ns=— OO 


(39) 


Since  the  incident  field  is  independent  of  angle,  the  scattered  pressure  is  also  angularly  inde* 
pendent  which  is  seen  in  Eq.  39.  The  signal- to-noise  ratio  resulting  from  an  isotropic  field  is 
subsequently 


(1p„  fti,. 


t  T.T-. 

nss— OO 


(40) 


Cosine  Squared  Directivity 

An  example  of  a  physically  significant  directivity  is  one  which  can  be  represented  as  a  di¬ 
rectional  cos*  function.*  The  directivity  function  in  cylindrical  coordinates  is  thus  prescribed 
by 


a  -f  cos*^', 
a  -j-  7Cos*^', 


— 7r/2  <<!>'<  ?r/2 
ir/2  <<l>’  <  37r/2. 


(41) 


The  Fourier  coefficients  for  this  specific  field  are  oo  =  a  +  0/2,  a_2  =  02  —  )8/4, 
60  =  a  -h  7/2,  and  6-2  =  62  =  7/4.  Applying  these  coefficients  to  Eq.  36  yields 


=  1,12a  +  (,3  +  ^)/2].  (42) 

K/i 


11 


Applying  the  same  coefficients  to  Eq.  37  yields 


=»|2c  +  (/3  +  7)/21  f;  T„r„- 


+  '(^  +  7)/4  E  +  e-^‘*T„TU] 

n=:— oo 


oo  oo 

n=s  — oo  r=:— oo 

[1  2  1  1  £ 

(n  — r  — 2)  n  —  r  (n  —  r  +  2) 


(43) 


resulting  in  a  signal-to-noise  ratio  of 


=  Y  + 


)3  +  7 


^  [e"''»T,r„V2  + 


2(4a  4-  /?  +  7) 

/?  00  00 

^  ^  2;r;e-dn-r)^(_l)(n-r-l)/2 

1 


7r(4a  +  ^  +  7)  r==^oo 

1  2 


+ 


(n  — r  — 2)  n  — r  (n  — r  +  2)J 


^n-~r,odd  • 


(44) 


Note  that  the  first  term  of  the  right  hand  side  of  this  equation  is  just  the  signal-to-noise  ratio 
of  the  isotropic  field. 


NUMERICAL  SIMULATIONS 

Measured  values  of  a  surface  directivity  led  to  a,  0,  and  7  values  of  0.33/75.4,  1.0/75.4, 
and  33.0/75.4  respectively.®  Our  numerical  simulations  were  carried  out  tising  the  MATLAB 
program  cyLinf  jig.elas  which  can  be  found  in  Appendix  F.  These  simulations  are  for  ka  val¬ 
ues  of  5,  20,  and  40  respectively  and  kR  values  of  16fca  and  32ka  which  represent  experimental 
measurement  distances. 
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Before  performing  numerical  simulations  on  an  elastic  cylinder,  the  impedance  term  of 
Eq.  32  is  expressed  in  terms  of  the  non-dimensional  variables  a/h,  c/cp,  and  pj p,  so  that 


n  Hn{kR)  Jill  .  2t  C  Nn{^) 

H'M  ^Cp  Dr^iQ) 


where 


(45) 


N„iQ)  =  _  n^), 


(46) 


and 


Dn(n)  =  _  n2)  _  £i£p 

pea 


0.*  -  fi2(l  +  /3V)  +  l3^n^ 


(47) 


For  our  simulations,  the  non-dimensional  variables  were  given  values  of  a/h  =  100,  Cpje  =  3.5, 
and  p,/p  =  7.5 

Since  we  deal  mainly  in  target  strength,  we  have  plotted  101og[(lps  |^}/(|p^  |^)]  with  respect 
to  angle  for  infinite  rigid  cylinders  and  infinite  elastic  cylindrical  shells.  In 
Figs.  1,  2,  and  3,  all  solid  lines  indicate  rigid  results  while  all  the  dashed  lines  indicate  elastic 
results.  For  each  figure,  the  outer  pair  of  curves  is  the  fimetion  evaluated  at 
kR  =  16*ka,  while  the  inner  set  of  curves  is  an  evaluation  at  kR  =  32*  ka.  For  high  ka,  the  de¬ 
crease  in  horizontal  scattering  {<f>  =  — 7r/2  and  <f>  =  'ir/2)  becomes  more  pronounced  and  in  both 
the  rigid  and  elastic  cases,  the  foward  scattering  from  ^  =  0  is  greater  them  the  backscattering 
from  this  direction.  In  a  comparison  of  the  elastic  and  rigid  scattering,  the  foward  scattering  is 
more  dominant  in  the  elastic  case  for  small  ka.  In  each  figure  it  is  also  seen  that  the  difference 
between  the  A;a  case  and  the  32  ♦  ka  case  is  on  the  order  of  3  dB  as  would  be  expected  when 
comparing  cylindrical  scattering  at  distances  which  differ  by  a  factor  of  two. 
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Fig.  1.  101og[P(^)]  for  the  cosine  squared  directivity  function  with  a,  P, 
and  7  values  of  0.33/75.4,  1.0/75.4,  and  33.0/75.4  respectively. 

^  indicates  the  direction  of  propagation  (solid)  and  the  direction 
from  which  the  wave  comes  (dashed).  The  center  of  the  plot  indi¬ 
cates  the  axis  of  the  cyhnder. 


SURFACE 


•■••ifVo)' 


♦*  ’’K 


Fig.  2.  101og[(|p5(^)/(|p;y  1^)1  for  infinite  rigid  cylinders  and  elastic  cylin 
drical  shells  evaluated  at  ka  =  5.  The  outer  solid(dashed)  line 
corresponds  to  the  rigid(elastic)  result  for  kR  =  16  *  ka.  The  in 
ner  solid(dashed)  line  corresponds  to  the  rigid(eiastic)  result  for 
kR  —  32  *  ka.  The  center  of  the  plot  indicates  the  axis  of  the 
cylinder. 


15 


SURFACE 
=  0 


_  A  A 

Fig.  3.  lOlogKIpgl  )/(|Pjv|  )]  for  infinite  rigid  cylinders  and  elastic  cylin¬ 
drical  shells  evaluated  at  ka  =  20.  The  outer  solid(dashed)  line 
corresponds  to  the  rigid(elastic)  result  for  kR  =  16  *  A:a.  The  in¬ 
ner  solid(dashed)  line  corresponds  to  the  rigid(elastic)  result  for 
kR  =  32  ♦  ka.  The  center  of  the  plot  indicates  the  axis  of  the 
cylinder. 


— 7r/2 
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Fig.  4.  101og[(|pj|  )/(|Pj»|  )]  for  infinite  rigid  cylinders  and  elastic  cylin 
drical  shells  evaluated  at  ka  =  40.  The  outer  solid( dashed)  line 
corresponds  to  the  rigid(elastic)  result  for  kR  =  16  *  ka.  The  in 
ner  solid( dashed)  line  corresponds  to  the  rigid(elastic)  result  for 


kR  =  32  ♦  ka.  The  center  of  the  plot  indicates  the  axis  of  the 
cvlind 


RIGID  SPHERE 


The  incident  pressure,  Pinc,  in  spherical  coordinates  for  a  plame  wave  propagating  in  the 
positive  z  direction  is  given  by'*-® 

Pinc{R,0)  =  Pi  exp\ikR  cos  $] 

oo  , 

=  Pi  X^(2n  +  l}i^Pn{cos9)j„(kR), 

n—0 


where  P,  is  the  amplitude  of  the  wave,  and  (R,d)  indicates  the  location  of  the  observation  point 
with  6,  the  azimuthal  angle,  measured  from  the  positive  z  eixis.  Pn  is  the  n'*  order  Legendre 
polynomial  and  jn  is  the  n**  order  spherical  Bessel  function.  Since  the  wave  is  propagating  in 
the  z  direction,  there  is  no  circumferential  angular  dependence.  The  scattered  pressure,  p,,  due 
to  this  incident  wave  striking  a  rigid  sphere  of  radius  a  is 

p,(R,e)  =  -Pi  f'(2n  +  (49) 

where  hn  is  the  n**  order  spherical  Hankel  function,  and  the  primes  indicate  derivatives  with 
respect  to  the  arguments. 

For  a  general  result  we  consider  a  plane  wave  propagating  in  the  {$',  <i>')  direction  which  is 
viewed  from  the  (6,  <f>)  direction  where  <f>  and  (f>'  represent  circumferential  eingles.  Let  0  indicate 
the  angle  between  the  two  directions  given  above.  The  relation  between  0  and  the  other  angles 
is 


cos  0  =  cos  0  cos  0'  +  sin  0  sin  0’  cos(0  —  0*).  (50) 

The  angle  0  now  replaces  the  0  of  Eqs.  48  and  49  so  that  the  Legendre  polynomials  become 
P„(cos0)  rather  than  P„(cos^).  Applying  Eqs.  3  to  5  of  Appendix  A  to  P„(cos0)  changes  the 


general  incident  and  scattered  pressures  to 


Pinc{R,^,  ,4>')  =  Piexp[ikR(cos6  cos6'  +  sin^sin^'  cos(o  — 


(51) 


and 

oo  n 

Ps{R,d,<j),e',(f>')  =  -P,  ^  ^  e,nX,..,nPn{c<^s6)PI^{cosB')coslmi<p  -  4>')]  (52) 

n=0  m=0 

respectively.  Xm,n  is  defined  to  be  to  be 


Ln,m 


=  (2n  +  l) 


s 

(n  +  m)! 


where  is  defined  in  terms  of  spherical  bessel  functions  as 


(53) 


rpS  •nJn(fcQ) 

”  Kika) 


hn(kR), 


(54) 


in  the  same  manner  as  the  rigid  cylindrical  case.  Cm  indicates  the  Neumeum  function  defined  in 
Eq.  77  of  Appendix  A.  From  this  point  we  suppress  the  S  superscipt. 

Equations  51  and  52  are  the  incident  and  scattered  pressure  due  to  one  plane  wave.  As  in 
the  cylindrical  formulation,  the  principle  of  superposition  is  used  to  generalize  to  the  situation 
of  a  spatially  continuous  noise  field.  The  continuous  incident  and  scattered  pressure  become 


(5',  ^')exp[ifci2(cos^cos5'  +  sin  ^  sin  0'  cos(<^  —  <^'))]dn' 


(55) 


and 


m  OO  n 


(56) 


n=0  ms=0 
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where  the  solid  angle  integration  dO,'  is  given  by  dCl'  =  sm9'd9'd0'  and  P^{9\  0')  is  the  direc¬ 
tional  amplitude  for  the  wave.  The  average  incident  pressure  squared  at  any  point  is  indicated 

by 


{\pA^,4>)\‘^)=  [  r s\n9'd9'd4>'  f  r\in9"d9"d<i>"{PAe',4>')P:{9"A")) 

Jq  Jo  Jo  Jo 

■  exp[i  kR{  cos  9  cos  9'  +  sin  9  sin  9'  cos(  <i>  —  4)'))] 

•  exp{— i/i;/2(cos^cos^"  +  sin  ^  sin  cos(  —  (t>"))]. 

Similarly  the  average  scattered  pressure  squared  at  ciny  point  is  indicated  by 


■  jr  f2tr 


(57) 


r  f  \m9'd9'd<f>'  r  f  \in9'^d9"dr{PAB\4>')P:{9",<f>’')) 
Jo  Jo  Jo  Jo 

oo  n 

■EE  emXn,mPn{cos9)PI^{cos9')cos[Tn{4>  -  4)')] 


n=0  m=0 
oo  j 


51  ^kXj,kPj  (cos 9)PjAcos  9")  cos[k{4>  -  4>")]. 

j=0  k=0 


(58) 


Let  the  correlation  between  the  pressure  amplitudes  be  defined  as 


(59) 


where  t>{9',  4>\  indicates  the  directivity  function  associated  with  the  noise  field.  As  in 

the  cylindrical  case  we  assume  that  the  pressure  fields  from  different  directions  are  imcorrelated 
so  that  the  directivity  function  reduces  to 


V{e\  4>\  4>")  =  4>')Hx)  =  V{9\  4>')^ 


sin^' 


(60) 
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Applying  these  definitions  to  Eqs.  57  and  58  yields 


IT  «27r 


jj 


sme‘d0'd<i>'v{e\  <t>') 


■  exp[ikR{cos6cos0'  +  sin 0 sin cos( -  o'))] 

■  expl—ikR{cos 6 cos 6'  +  sin^sin^'  cos(<j>  — 


nsm6’de'd<i>'V{e',(f>'), 


(Cl) 


and 


IP.V 


oo  n  oo  i 

=  EEEE  €m€ltXn.mX2.fc-Pr(cose)P*(cos^)/(n,m,j,  Ir), 

n=0  mssO  j—O  k=Q 


(62) 


where 


I{n,m,j,k)  = 


d6'd<i>'  sin  ^'P^’Ccos  ^')P/(cos  e*)V{0\  <p') 


•  cos[m(^  —  <^')]cos[t(^  —  4>')]. 


(63) 


CIRCUMFERENTIALLY  SYMMETRIC  DIRECTIVITY 

In  the  case  that  the  directivity  is  circumferentially  symmetric,  V{9\<f>')  —  T>(0'),  the  result¬ 
ing  average  noise  field  is 

=  27r  /  V(0')sm0'd0'.  (64) 

I-*  7  1  Jo 

After  the  simplifications  of  Appendix  C  have  been  made,  the  scattered  pressure  field  reduces  to 


IP,  I' 


OO  OO 

=  2ir^^(2n  +  l)(2j  +  l)T„r;  E 

n=0  j=0  m—0 


(65) 
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where 


Sn,j,m  — 


(n  -  m)!  (>  -  m)! 


(„  +  ;^y.  (  .  j.n,: 


(66) 


and 


sin  d'P^^’Ccos  0')P/(cos  e')Vie')de'. 


(67) 


For  every  pair  (nj),  there  is  also  a  pair  {j,n)  which  gives  the  same  Sn,},m-  However,  TnTJ  be¬ 
comes  TjT*.  By  introducing  a  change  to  the  summation  limits,  the  expression  for  the  scattered 
pressure  becomes 


(Ip.  I') 
\PA' 


"■E  +  mi  +  i)ir„r;  + 

nstO  ji=0  m=0 


(68) 


Cosine  Squared  Directivity 

We  define  V(ff')  as  we  did  in  the  case  of  the  cylinder,  but  give  it  circumferential  symmetry. 
Therefore 


a  +  j3  cos^  $\ 
a  +  7cos^0', 


0  <  0'  <  7r/2 
7r/2  <  B'  <  TT. 


(69) 


The  average  noise  field  for  this  directivity  becomes 


=  27r 


/•»  /■» 
a  I  sin  6'dB'  +  y?  /  cos^  B'  sin  B'dB'  +  7  /  cos^  B'  sin  B'dB' 
Jo  Jo  Jit/2 


=  2?r[2a  +  (/9  +  7)/3]. 


(70) 
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After  the  simplifications  of  Appendix  E  have  been  made,  the  average  scattered  pressure  be¬ 
comes 


=  4^a  '£(2n  +  1)T„T„- 
IP,  I  ^ 


n=0 


n=0  m=0 


+  2x{^  +  7)  •  s 


n=2 


n-2 


m=0 
n  — 1 


n=l 
/  00  n  — 1 


m=0 


EE  44TvT;+TjT:]  Yl  emP„""‘(cos(9)P-”*(cos0) 

n=l  jsaO 


m=0 


(n  -  m  +  l)(i  -  m  +  l)P„+ij+i,m 


00  n  — 1 


J-1 


+  2jr(;3  -  7)  •  S 


+ E  E  + p>p»'i  E  ‘».p„-"(cos«)Pj-'"(cos«) 

n=2  j  =  l  »n=0 

•  (n  -  m  +  1)0'  +  m)Pn+ij-i,m 


00  n— 1 


+  ^  wlTnr;  +  r,r;|  53  e„F,-"(cos(i)p-”(cos«) 

n=l j=0  m=0 

•  (n  +  m)(j  -  m  +  l)P„-i,>-n,m 

+  E  E  +  Ti-n]  «„p-”(cos«)P-’"(cos») 

n=2  jssl  m=0 

.  (n  +  m){j  +  m)Vn-i,j-i,m- 


(71) 
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Dividing  Eq.  71  by  Eq.  70  yields  the  proper  expression  for  the  signai-to-noise  ratio  for  a  cir¬ 
cumferentially  symmetric  cosine  squared  directivity. 

NUMERICAL  SIMULATIONS 

As  with  the  cylinders,  we  have  plotted  101og[(lp5l^)/(|p^  j^)]  with  respect  to  azimuthal  an¬ 
gle  for  a  rigid  sphere  in  a  circumferentially  symmetric  noise  field.  We  use  the  samecoefficients 
of  the  cosine  squared  directivity  so  that  a  =  0.33/75.4,  ^  =  1.0/75.4,  and  7  =  33.0/75.4.  it 
should  be  noted  that  these  values  of  a,  0,  and  7  result  in  a  normalized  noise  condition  which 
reduces  Eq.  70  to  unity.  Our  numerical  simulations  were  carried  out  using  the  MATLAB  pro¬ 
gram  sphjig-dir  which  is  in  Appendix  F.  The  simulations  are  for  ka  values  of  5,  20,  and  40  re¬ 
spectively  and  kR  values  of  16A:a  and  32ka. 

In  Figs.  4,  5,  and  6,  all  solid  lines  indicate  the  function  evaluated  at  kR  =  16*  ka  while 
dashed  lines  indicate  an  evaluation  of  the  function  at  kR  =  32  *  ka.  For  high  ka,  the  decrease 
in  horizontal  scattering  (6  =  7r/2)becomes  more  pronounced  and  the  foward  scattering  from 
^  =  0  is  greater  than  the  backscattering  from  the  same  direction.  In  each  figure  it  is  also  seen 
that  the  difference  between  the  16*  ka  case  and  the  32*  ka  case  is  on  the  order  of  6  dB  as  would 
be  expected  when  comparing  spherical  scattering  at  distances  which  differ  by  a  factor  of  two. 

CONCLUSION 

An  expression  was  obtjiined  for  the  signal-to-noise  ratio  due  to  scattering  from  an  infinite 
rigid  cylinder  or  an  infinite  elastic  cylindrical  shell  subjected  to  a  normally  incident  ambient 
noise  field.  This  solution  is  applied  to  the  specific  case  of  isotropic  fields  and  to  the  case  of 
noise  fields  which  can  be  expressed  by  a  squared  cosine  function.  Similar  methods  were  also 
applied  to  obtain  the  signal-to-noise  ratio  due  to  the  scattering  from  a  rigid  sphere  subjected  to 
such  a  noise  field  with  the  assumption  of  circumferential  symmetry. 

The  effect  of  the  ambient  field  is  to  smooth  the  scattering  function  which  results  from  a 
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single  incident  plane  wave.  There  is  still  a  larger  foward  scattering  than  back  scattering,  but 
the  typical  lobes  observed  from  one  incident  plane  wave  are  no  longer  distinguishable  due  to 
the  contributions  from  all  the  directions. 

The  limitations  of  these  results  arise  from  the  assumption  that  the  noise  field  in  the  cylin¬ 
drical  case  is  assumed  to  be  normally  mcident  to  the  cylindrical  surface.  It  is  expected  that 
similar  procedures  can  be  applied  to  obtain  results  for  a  more  realistic  situation.  It  is  also  ex¬ 
pected  that  similar  mathematical  techniques  can  be  apphed  to  other  directivity  functions  and 
to  the  case  of  an  th\stic  spherical  shell. 
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Fig.  6.  101og[(|pjp)/(|pjy  1^)]  for  a  rigid  sphere  evaluated  at  ka  =  20. 

The  solid  line  corresponds  to  kR  =  16  *  it;a  and  the  dashed  line 
corresponds  to  kR  =  32*  ka. 
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APPENDIX  A 


MATHEMATICAL  IDENTITIES 


INTEGRALS  OF 


d4>'  =  27r6,,o 


(72) 


(73) 


and 


/•3t/2  9 

/  d<i>'  =  7r^,,o  -  . 

Jir/2  V 


(74) 


PROPERTIES  OF  ASSOCIATED  LEGENDRE  FUNCTIONS^ 


Pn{cos  V>)  =  Pn(cos  6  COS  +  sin  6  sin  $'  cos(<^  —  <!>*)) 


=  S  «m7^^37^^r(cos^)p;^(cos0')cos(m(^-^')|, 
(n  +  m): 

«nssA  ^  ' 


(75) 


where  are  the  associated  Legendre  functions  defined  by 


Pr(x)  =  (-1)-(1 


(76) 
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and  €m  indicates  the  Neumann  function  defined  by 


f  1,  if  m=0; 

\  2,  otherwise. 


(TT) 


(78) 


n  ,7  5 


(79) 


=  ,,  ^  m  +  l)PZ-i(.=‘)  +  (n  + 

(2n  +  Ij 


(80) 


and 


i-ir 


(n  —  m)! 
(n  +  m)! 


(81) 
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APPENDIX  B 


EVALUATION  OF  THE  SCATTERED  PRESSURE  DUE  TO  A  CYLINDER 
IN  AN  AMBIENT  FIELD  WITH  A  GENERAL  DIRECTIONAL  DIRECTIMTY 


The  integrations  of  Eqs.  34  and  35  can  be  evaluated  using  Eqs.  73  and  74  of  Appendix  A. 
Applying  these  to  Eq.  34  with  rj  =  m  results  in 

{\Ps\")  ^  ^  .  f_c  , 


Om  I  ^  dm, odd] 

|■‘/l  m=-(X)  ^  ^ 

"I"  ^  ^  dm  2 


'm,odd 


m=  —  oo 


~  (_l)(m-l)/2 

=  jr(ao  +  fco)  +  2  2__^  —  dm) - 3;; - dm, odd 


(82) 


m 


Applying  the  same  equations  to  Eq.  35  with  tj  —  n  —  r  +  m  yields 
{\Ps\") 


00  00 


\PA 


=  E  E 


— i(n— r)^ 


na=— 00  — 00 

<x> 


^  /  (_l)(n-r+m-l)/2  \ 

/  ^vn  I  ^*T»— r+m,0  4"  2  ■  . .  r+m,o({ii  1 

\  n  —  r  +  m  J 

l=— oo  '  ' 


Lm=— oo 
oo 


r  +  m 

—  r+m  — 1)/2 


4"  /  .  dm  I  ^^T»— r+m,0  2  '  dn  —  r+m,odd  1 

'  V  n  —  r  +  m  } 

m=**-oo  '  - 


Using  the  fact  that  6n-r+m,o  =  dn-r,~mi  Eq-  83  simplifies  to 


{\Ps\") 

\PA" 


=7r  +  J]  r„r; 


n+m 


m=*oo 

oo 


n=s— oo 
oo  oo 


‘r 


+  2  ^  (a„-6„)  53  T.T, 

OO  n=  — oors^-oo 

g-t(n-r)0^_2^(n-r+m-l)/2 
(n  —  r  +  m) 


dn— r+m, odd  • 


(83) 


(84) 
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APPENDIX  C 

EVALUATION  OF  THE  SCATTERED  PRESSURE  FROM  A  RIGID  SPHERE 
DUE  TO  A  GENERAL  CIRCUMFERENTIALLY  SYMMETRIC  DIRECTIVITY 


As  a  result  of  circumferential  symmetry,  I  can  be  sepai'ated  in  the  following  manner. 


where 


h'{n,m,j,k)=  [  sme'P;;^{cos6')Pj^(cos0')V{e')d9' 
Jo 


,2» 

k)  =  I  cos[m{<f>  —  (f>')]  cos[A:(i^  —  J>')]d<f>' 

Jo 


Rewriting  the  integrand  of  k)  as 


cos[m(^  —  4>’)\  cos[fc(0  —  ^  e'^^e  , 


where 


^  =  {m  +  A:,  — {m  +  Ar),m  —  fc,  — (m  —  A:)} 


and  applying  Eq.  72  of  Appendix  A  yields 


/^(m,  fc)  =  j  ^ 

TT 

— (m  — fc),o) 

“  ‘^(^m,0^k,0  "F  ^m,k)' 
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Therefore,  for  a  circumferentially  symmetric  noise  field,  the  scattered  field  from  a  rigid  sphere  is 
obtained  by  applying  this  separation  to  Eqs.  62  and  63.  The  result  is 


oo  oo 

-EE  Xn,oXj,oPnicos9)Pjicos  9)Ie-{n,  0,  j,  0) 

n=0  j~0 


oo  oo  m«n[n,j] 

+'EE  E  ^mXn,mXj,mPn(cos9)Pp{cos9)Ie'in,m,j,m} 

n=0  7=0  m=0 

oo  oo  min{n,j] 

=  2n  EE  E  «mX:n.mXlm-Pr(cos6)P,^(cose)7fl-(n,m,;,m). 

n=0  7=0  m=0 


(91) 


Rewriting  this  in  terms  of  T„ 


|e,l" 


oo  oo  mmfnj) 

2ir'^'£(2n  +  l)(2i  +  l)T„T;  E 

n=0  7=0  m=0 


(92) 


where 


'n,7,m 


(n  —  m)!  (j  —  m)! 
(n  +  m)!  (j  +  m)! 


P^(cos  9)Pp(cos  9)Ig'{n,  m,j,  m). 


(93) 
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APPENDIX  D 


INTEGRATION  OF  TWO  ASSOCIATED  LEGENDRE  FUNCTIONS 
OVER  AN  ARBITRARY  RANGE  CONTAINED  IN  [-1,1] 


Consider  the  integral 


P-(x)Pf(i)dx 


2  (n  +  m)! 
2n  +  1  (n  —  m}! 


(94) 


where  m  <  n  and  m  <  j.  This  is  a  well  known  integral,  but  when  the  integration  limits  are 
different,  the  integration  is  not  zis  easily  determined.  Therefore,  we  have  derived  an  expression 
for  general  integration  limits. 

The  Associated  Legendre  fimctions  satisfy  the  differential  equation® 

(1  -  x^)t/"  -  2xy'  +  |n(n  -f- 1)  -  ^  |  y  =  0 


which  is  equivalent  to 

~{(l-x2)y'}  +  |n(n  +  l)-Y-^^|y  =  0.  (96) 

Multiplying  by  (1  —  x^)  and  substituting  P^  and  Pp  for  y,  Eq.  96  becomes 

(1  -  {(1  -  x^)Pr'}  +  {(1  -  x")n(n  +  1)  -  m)  Pr  =  0  (97) 

and 

(1  -  {(1  -  x")P”' }  +  ((!-  +  1)  -  m)  P”  =  0,  (98) 
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where  the  primes  indicate  derivatives  with  respect  to  x  and  the  argument,  x,  of 
suppressed.  Multiplying  Eqs.  97  and  9S  by  PJ"  and  P^  respectively  yields 

(1  -  {(1  -  }  P/"  +  1(1  -  +  1)  -  m)  =  0 

and 

(1  -  x2)A  {(1  _  x^)pp']  p;r  +  {(1  -  + 1)  -  m)  prPn  =  0. 

When  j  aLid  n  axe  different,  we  subtract  Eq.  100  from  Eq.  99  and  factor  the  (1  — 
obtmn 

^  {(1  -  I*)P„"' }  P”  -  £  {(1  -  }  Pn  =  U(j  +  1)  -  +  1)) 

An  application  of  the  chain  rule  to  the  quantities  (1  —  yields 


^  {(1  -  x^)pfpr} = p”^  {(1  -  x^)pr'} + (1  -  x')p„"’pr' 


and 


£  {(1  -  x*)P”Pf '}  =  p„”^  {(1  -  x^ipf' } + (1  -  x")p"'pr'' 


Upon  subtracting  these  equations  we  obtain 


dx 


{(1  -  T^)pz‘pr}  -  s  {(1  -  =‘^)p"pr’ }  = 


and  Pf  is 

(99) 

(100) 

)  term  to 

'■  (101) 

(102) 

(103) 

(104) 
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Substituting  Eq.  104  into  Eq.  101  yields 


■£  {{1  -  x^)[Pfpp'  -  PTPf\]  =  W  +  1)  -  n(r.  +  1))  PfPr- 


(105) 


Integrating  this  equation  over  the  interval  [a,  6]  which  is  contained  in  the  interval  [—1.1]  yields 


f  PZ(x)P^{x)dx 

J  a 


[P^\x)P^(x)  -  P„”-(x)P,r"'(x)]|; 
j(j  +  l)-n(n  +  l) 


(106) 


as  a  general  expression  for  the  integral  of  the  product  of  associated  Legendre  functions.  When 
a  =  0  and  6=1  this  reduces  to 


[p-(0)Pf'(0)-pr'(0)^r(Q)] 

j(j  +  1)  -  n{n+  1) 


(107) 


It  should  be  noted  that  this  integration  is  true  for  m  <  n  ^md  m  <  j.  If  m  >  n  or  m  >  j,  then 
the  integral  is  zero.  It  should  also  be  noted  that  when  n  =  j,  Eqs.  99  and  100  are  the  same.  If 
they  were  subtracted  to  obtain  Eq.  101,  the  resulting  equation  would  be  0  =  0.  Therefore  when 
n  =  j  the  integral  must  be  evaluated  in  a  different  manner. 

The  integrand  of  Eq.  94  is  even  so  that 


P-(z)Pr(x)dx 


1  (n  +  m)! 
2n  +  1  (n  —  m)! 


(108) 


Also,  for  n  ^  j  Eq.  94  can  be  rewritten  as 

j'  P^{x)P^{x)dx  =  [1  +  (-1)"+^]  p;p{x)pp{x)dx  =  0. 


(109) 


Therefore  when  n  +  j  is  an  even  number 

/’  p”(i)P7(iMx  =  0. 

Jo 


(110) 
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Using  these  facts,  Eq.  107  can  be  rewritten  more  precisely  as 


I 


(111) 


where 


n,j,m 


P^(0)Pf'(0)  -  P^‘(0)P,’"(0} 

j(j  +  l)-n(n  +  l) 


and 


f  1,  0  <m  <  M; 
^  0,  otherwise. 


(112) 


(113) 
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APPENDIX  E 


EVALUATION  OF  THE  SCATTERED  PRESSURE  FROM  A  RIGID  SPHERE  DUE 
TO  A  CIRCUMFERENTIALLY  SYMMETRIC  SURFACE  NOISE  DIRECTIVITY 


(114) 


Applying  Eq.  69  to  Eq.  67  reduces  the  integral  term  for  the  scattered  pressure  to 

Ie'(n,  m,j,m)  =  a.  f  sm(6')PJP{cos6')Pp{cos6')d6' 

Jo 

fitl2 

+  0  sin(^')cos2(e')P^(cos0')JT(cos^')d0' 

Jo 

+  'r  f  sin(^')cos2(^')P;^(cos5')^r(cos^')d^' 

Jr{2 

=  a  £  P™(x)P”(xMx  +  [^  +  (-1)"+M  [  x^PZMP"{x)di 


The  latter  expression  is  obtained  using  the  substitutions  x  =  cos(^'),  dx  =  —  sin{6')d6' ,  and  ap¬ 
plying  Eq.  78  of  Appendix  A.  For  the  first  integration,  apply  the  orthogonality  relation  of  Asso¬ 
ciated  Legendre  ftmctions  as  given  in  Eq.  79  of  Appendix  A,  and  in  the  second  integral,  apply 
the  identity  given  in  Eq.  80  to  xP^(x}  and  xPp(x).  The  second  integral  therefore  becomes 


f'  = 

Jo 


(2n  +  l)(2j  -f  1) 


(n  -  m  -f  l)(j  -  m  -I- 1)  f  Pn^i{x)PJli{x)dx 

Jo 

-I-  (n  -  m  -I-  l){i  +  m)  /  P,^i(x)P/!li(x)dx 

Jo 

-f  (n-}-m)(j -m-f-1)  /  P„"Li(a;)P;!;i(x)dx 

Jo 

+  {n  +  m){j  +  m)  f  Pn-i{^)Pjliix)dx 

Jo 


(115) 
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In  Appendix  D  it  was  shown  that, 


I 


P™(x)Pf (xM^:  = 


(116) 


When  this  is  applied  to  the  integrations  of  Eq.  115,  the  second  integral  of  Eq.  114  becomes 


f  x^p;;^{x)pp{x)dx  = 

Jo 


(2n  +  l)(2i  +  l) 


(  fn  +  l)(j  ^  "b  l)^n+l,j-t-l,m<^n4-J  +  2,od</^m,n+l bnij  +  l  ^ 

+  (n  -  m  +  l)(i  +  m)'P„+ij_j  ddUm 

"b  ip  "b  ^)(j  m  +  l)^n— l,m^n+j,o  ddU  m,n— 1  UmJ+i 

I  "b  "b  +  J72)^n— 1  ,j  — l,m<5n+j  — 2,o<l<i^m,n— 1  ^m,)  — 1  •  J 


(117) 


+ 


(n  -  m  +  l)(i  —  m  +  1)  (n  +  m  4- 1)! 


(2n  4-  3)  (n  -  m  + 

(n  -  m  4-  l)(j  4-  m)  (n  +  m  4- 1)! 


(2n  4-  l)(2i  4- 1) 


4- 


4- 


+ 


(2n4-3)  (n  —  m  4-1)! 

(n  4-  rn)(j  —  m  4- 1)  (n  4-  m  —  1)! 
(2n  —  1)  (n—m  —  1)! 

(n  4-  m)(j  4-  m)  (n  4-  m  -  1)! 


U„, 

,Tl+l 


(2n  —  1)  (n  —  m  —  1)! 


*^n— l,ji  — — 1 


Applying  this  and  the  delta  function  identities  6n±i  j±i  =  =  ^n±2, ji 
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^n+jiL2,odd  —  ^n+j,odd  reduces  Iff’  tO 


r  .  •  ^  2  (n  +  m)!. 


(  (i  -  +  1)  (n  -f  m  -f  1)! 

(2n  +  3)  (n  —  m)! 


+ 


(;3  +  7] 


(2n  +  l)(2i  +  l) 


^n,j^m,n+i 


.  (j  +  m)  (n  +  m-M)!^  rr 

+  TT - ^ ; - r: — o„4.2  iU, 


+ 


r2n  +  3)  (n-m)l 
(J  -  m  +  1)  (n  +  m)! 


(2n  —  1)  (n  —  m  —  1)! 


^n— 2,j  — 1 


,  0  +”j)  (n  +  m)!  ^ 

T"  71^ - 7T7  TTr^n,j^m,n  —  J 

(2n  —  1)  (n  -  m  —  1): 


~  7]^n+j,o(f<< 
(2n  +  l)(2j  +  l) 


'  (n  -  m  +  l)(j  -  m  +  l)'Pn+l,>  +  l,mt^m.n  +  ll^m,>+I 
+  (n  -  m  +  1)0’  +  Tn)V„+ij-i,ml^m,n+l^mJ-I 

1 

+  (n  +  m){j  —  m  +  l)^n— l,j+l,mt^m,n— 

.  +  (n  +  Tn){j  +  — — 1 


> 


(IIS) 
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Substituting  this  expression  into  Eqs.  66  and  68  yields 


\P.t 


=  47ra  'Y^{2n  -I-  l)r„r;  c, 


nssO 


m=sO 


(n  +  m)\ 


+  '1/3  +  7]  53  E  H.-dTni;  +  T,TX\ 

n=0  j=0 

(n  —  m)!  (j  —  m)! 


E  I  i  ri”  iOcos^)P-(cos^) 

(n  +  m)!  0  +  m)!  " 

m=0  '  ' 


(  (j  -  m  4- 1)  («  +  m  -f  1)! 
(2n  +  3)  (n-m)! 


,  (j+m)  (n  +  m  +  1)!^ 

(2n  +  3)  (n-m)! 

.  O’-m  +  l)  (n  +  m)!  ^ 

+  (2n-l)  (n-m- if 

.  0+m)  (n  +  m)!  ,  „ 

I  (2n  -  1)  (n  -  m  - 


oo  n 

+  -  T]  ^  6n+i,od<lHn-j)[TnT-  +  TjT*] 

n=0  7=0 


E  i"  ""i!  ifr^v  g)pr  (cos ») 

^  (n  4  m)!  {j  4  m)\ 

m=0  '  /  v-” 

(n  —  ni  4  l)(j  >Tl  4  l)^n4-l,j+l,Tn^n»,n+l^ni,y4'l  ^ 
4  (n  —  T7Z  4  1)0"  d"  —  —  \ 


4  (n  4  m)(j  —  m  4  l)7^n-i  j+i.m  Um,n-lU,n  ,>+l 
.  4  (n  4  fn){j  +  ^)'Pn-l,j-l,mVm,n—lUTn,j-i 


(119) 
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Note  that  whenever  n  +  j  is  odd,  e„_j  =  2,  (  — I)"'*’-'  =  —1  and  the  limits  n  =  [0, ooj.  j  =  (0. n] 
can  be  replaced  by  n  =  [l,oo],  j  =  [0,  n  —  1).  Using  these  identities,  applying  Eq.  81  of  Ap¬ 
pendix  A  to  the  ±  7  terms,  and  applying  Eq.  75  of  Appendix  A  to  the  q  term  simplifies 
Eq.  119  to 

=  47,0  f^{2n  +  1)T.T„- 
\“l  I  n=0 


oo  n  3 

+  »l/9  +  7l  +  Y. 

n=0j=0  m=0 

f  (j  -  m  4- 1)  (n  H-m -f  1)!  c  tt 
(2n  +  3)  (n-m)! 

.  (j  +  m)  (n  +  m  +  1)!  ^ 

+  /«  I  o\  ■  -  On+2,>‘''m,n-(-l 

(2n  +  3)  (n  —  m)! 

.  <  . 

.  (i-m  +  1)  (n  +  m)!  ^ 

(2n-l) 


V 


+ 


(i+m)  (n  +  m)!  ^  „ 

'}n' . ■"7\  /  1 

(2n  —  1)  (n  —  m  —  1)! 


/ 


(120) 


+  27,1/3  -7lYY  +  ^7^n]  E  "(cosdjP-^CcosS) 

n=l  j=0  ms=0 

'  (n  -  m  +  l)(j  -  m  +  l)P„+ij+i,mi^m,n+lUm,j+l  ' 

+  (n  -  m  +  l)(ji  +  m)P„+ij-i^rn^m,n+ll^mJ-l 
+  (n  +  Tn)(j  -  m  +  l)P„_i,>+i,mt^m,n-lt^m,j+l 

,  +(n  +  m)(j +m)P„-ij„i,n,Um,n-l^m,3-l 
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Eliminating  the  delta  functions  and  step  functions  yields  an  average  scattered  pressure  of 


il||^  =  4^0  f;(2n  +  l)r„r* 


n=0 
(  oo 


n=sO  m=0 

oo 


+  521T„T„-_2  +  T„.jT„*l 


+  27r(^  +  7)  •  S 


n~2 


m=0 
n  — 1 


2)! 


+ Ei^nr;]  E 


n=l 
(  00  n— 1 


m=0 


Y,  E ^"+j.-'4r„T;  +  Tjn\  E  «™e-”(cos«)p-'"(cos») 

n=l  jssO 


msO 


(n  -  m  +  l)(j  -  m  +  l)Pn+i,>+i,m 


00  n  — 1 


+  EE  5n+i,o  MTnT;+TiT:]  E  e„P„-"(cos«)p-'"(cos«) 

n=2  j=l 


+  27r(/3  -  7)  •  S 


m=0 

•  (n  -  m  +  l)(i  +  m)7^„+i  J_i,r 

00  n— 1  j 


+  E  E  +  I>T„’l  E  <„P„-'"(cos«)P,-"'(cos9) 

u=l j=0  m=0 

•  (n  +  m)(j  -  m  +  iyPn-i,j+i,m 

+EE  ^n+>,o  +  Tir;i  E  t»p,-"(co8e)p-”‘(cos9) 

n=s2  j=l  n>=l) 

.  (n  +  m)(j  +  myPn-  l,m« 


(121) 
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APPENDIX  F 

MATLAB  PROGRAMS  FOR  EVALUATING  THE  AMBIENT 
NOISE  SCATTERING  FROM  INFINITE  RIGID  AND  ELASTIC 
CYLINDERS  AND  FROM  RIGID  SPHERES 
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oVO  0>P  o\P  o'? 


cyl_inf_rig_elas 

This  program  determines  the  scattering  from  an  infinite  rigid 
cylinderand  an  infinite  elastic  cylinder  due  to  a  directional 
%  ambient  noise  field  described  by  the  directivity 
%  D(PHI)  =  Sum[a(m)  e:>q:)(iin*phi)  ] ;  -pi/2  <  PHI  <  pi/2 

%  D(PHI)  =  Sum[b(m)  e>^(im*phi)  ] ;  pi/2  <  PHI  <  3pi/2. 

%  The  above  directivity  function  is  a  general  fourier  ejq^ansion, 
%  and  the  program  is  written  so  that  any  directivity  of  this 
%  form  can  be  used.  Specifically 

%  D(PHI)  =  alpha  +  beta*cos (PHI) ^2;  -pi/2  <  PHI  <  pi/2 

%  D(PHI)  =  alpha  +  garTma*cos  (PHI)  ^2;  pi/2  <  PHI  <  3pi/2 

%  is  used.  Frcm  this  directivity  fimction  a  (m)  and  b  (m)  are 
%  a(0)  =  alpha  +  beta/2  a  (2)  =  a  (-2)  =  beta/4 

%  b(0)  =  alpha  +  gamma/2  b(2)  =  b(-2)  =  gaiTma/4. 

% 

% 

% 

clear 

% - FOURIER  COEFFICIENTS - 

alpha  =  0.33/75.4; 
beta  =  1.0/75.4; 
gamma  =  33.0/75.4; 

mvec  =  [-2  0  2];  %  it  for  nonzero  coeff. 

am  =  [beta/4  (alpha  +  beta/2)  beta/4];  %  [a (-2)  a(0)  a (2)] 

bm  =  (gamma/4  (alpha  +  ganina/2)  ganiTia/4];  %  [b(-2)  b(0)  b(2)j 
=  am  +  tm; 
airb  =  am  -  tm; 

% - ^WAVENUMBER  VARIABLES - 

ka  =  20.0;  %  nondimensional  wavenuirber 

kamult  =  [16  32] ;  %  observation  points  are  at 

kR  =  ka*kamult;  %  one  and  two  boat  lengths 

kRka  =  [kR(l)  kR(2)  ka]; 
numka  =  length  (kRka) ; 


%  - msE  vmjES - 

phimin  =  -pi/2.0; 
phiinax  =  3*pi/2; 
nurrphi  =  201; 

phistep  =  (phimax-phiitiin)  /  (nurrphi-1) ; 

phi  =  phimin  :phistep:phimax;  %  angle  vector 

%  - ELASTICITY  CONSTANTS - 

ccp  =  1/3.5;  %  c/cp  —  sound  speed  ratio 

rhosrho  =  7.500;  %  rhos/rho  —  density  ratio 

ha  =  1/100;  %  h/a  —  thickness/radius  ratio 

Omega  =  ka*ccp;  %  nondimensional  freq.  parameter 

betasqr  ==  (ha'^2)/12; 

% - PPOGRAM  PREPARATION - 

maxnr  -  30;  %  sum  to  maxnr  rather  than  inf. 

maxm  =  max  (abs  (mvec) ) ; 
iplusn(0+l)  =  1; 
for  n  =  l:maxnr+maxm 

iplusn(n+l)  =  i*iplusn(n); 
end  %for  n — 1 
iirphi  =  i*invec. '*phi; 
eiiiphi  =  ej^(ir[phi) ; 

% - ^BESSEL  CALL - 

[ J/  dJ,  y,  dY]  =  vJdJYdY  (maxnr+maxm,  kRka) ;  %  J  (n,  k)  =J  (n,  kR  (k) ) ; 

H  ==  J  +  i*Y;  %  J(n,numka)  =  J(n,ka) 

dHka  =  dJ(:,nuraka)  +  i*dY { : , numka) ;  %  only  need  H’  (ka) 

%  - ^ELASTICITY  TERMS - 

for  n  =  Oimaxnr-fmaxm 
nsquare(n+l)  =  n''2; 
nfourth(n+l)  =  n''4; 
nsixth(n+l)  =  n''6; 
end  %for  n — 2 
cm2n2  -  Qraega‘'2  -  nsquare; 
num  =  cin2n2*2*i*ccp/pi; 
deni  =  Qnega*H(:, numka)  .*an2n2.'; 
den2  =  ha*rhosrho*dHka/ccp; 


den3  =  Cn^a'"4-Qmega''2*  (l+nsquare+betasqr*nfou2:th)  +betasqr*nsixth; 
den  =  denl-den2.*den3. ' ; 
elastic  =  num. ' . /den; 

% - PROGRAM - 

for  R=l:numka-1 

coef  f n_r  =  iplusn . ' .  *dJ  ( : ,  nurnka) .  *H  ( :  /  R) .  / dHka; 
coeffr_r  =  conj (coeffn_r) ; 

coeffn_e  =  iplusn. ' . * (dJ ( : , numka) +elastic) . *H ( : , R) . / dHka; 
coeffr_e  =  conj (coeffn_e) ; 

sumlin_r  =  zeros  (phi) ; 
sumlm_e  =  zeros  (phi)  ; 
for  j  =  1: length (mvec) 
m  =  invec(  j) ; 
sumn_r  =  0; 
sumn_e  =  0; 
for  n  =  ^naxnrrmaxnr 

suinn_r  =  sunin_r  +  coeffn_r  (abs (n)  +1)  *coeffr_r (abs  (n+m)  +1) ; 
surnn_e  =  sumn_e  +  coef fn_e  (abs  (n)+l)*coeffr_e  (abs  (n+m) +1 ) ; 
end  %for  n — 3 

sumln^r  =  sumlit^r  +  apb(j)  *eiiiphi  (j, :)  *sunin_r; 
sumlm_e  =  sumlm_e  +  apb(j)  *eiirphi(j/ :)*sumn_e; 
end  %for  j — 1 
sumln^r  =  pi*sumln^r; 
sumlin_e  =  pi*sumlm_e; 

sum2m_r  =  zeros  (phi) ; 
sum2m__e  =  zeros  (phi) ; 
for  j  =  1:  length  (mvec) 
m  =  mvec(j) ; 
surnnr_r  ==  zeros  (phi) ; 
suinnr_e  =  zeros  (phi); 
for  n  =  “inaxnr:maxnr 
for  r  =  -maxnrimaxnr 
odd  =  (-1)  (n-r4m) ; 
if  odd  =  -1 

enrphi  '  ejp  (-i*  (n-r)  *phi) ; 

cMnr_r  =  coef  fn_r  (abs  (n)  +1)  *coeffr_r  (abs  (r)  +1)  ; 
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sumnr_r  =  sumnr_r  +  (-1) ( (n-r+m-l) /2)  *chinr_r*enrphi/ (n-r+m)  ; 
chinr__e  =  coeffn_e  (abs  (n)  +l)  *coeffr_e  (abs  (r)  +1) ; 
surnm:_e  =  suiTinr_e  +  (-1)  ^  ( (n-r+rrt-1) /2)  *chinr__e*enrphi/ (n-r+m)  ; 
end  %  if — odd 
end  %for  r — 1 
end  %for  n — 4 

suni2m_r  =  sum2m_r  +  aiTto(  j)  *surnnr_r; 
sum2m_e  =  suin2m_e  +  aiTb(  j)  *sumnr_e; 
end  %for  j — 2 
sum2m_r  =  2*sum2m_r; 
sum2m_e  =  2*sum2m_e; 

sigtonoise_r  (R, : )  =  (surnlin_r  -f  sum2m_r)  /  {pi*af±)  (2) ) ; 
sigtonoise_e  (R, : )  =  (sumlm_e  +  sum2m_e) / (pi*apb (2) ) ; 

TS_r  (R, : )  =  10*logl0 (sigtonoise_r  (R, : ) ) ; 

TS_e  (R, : )  =  10*logl0 (sigtonoise_e (R, : ) ) ; 
end  %for  R — 1 
end 
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%  This  program  determines  the  scattering  fron  a  rigid  sphere  due  to  a 
%  directional  airfoient  noise  field  described  by  the  directivity 
%  D (THETA)  =  alpha  +  beta*cos (THETA) ^2;  0  <  THETA  <  pi/2 

%  D  (THETA)  =  alpha  +  gamma*cos  (THETA)  ^2;  pi/2  <  THETA  <  pi 

% 

% 

clear 

%  - ^DIPECTIVITy  COEFFICIENTS - 

alpha  =  0.33/75.4; 
beta  =  1.0/75.4; 
gamma  =  33.0/75.4; 

% - WAVENUMBER  VARIABLES 

ka  =  5.0; 
kamult  =  [16  32]; 
kR  =  ka*kamLLlt; 
kRka  =  (kR(l)  kR(2)  ka]; 

% - WANGLE  VALUES - 

numtheta  =  101; 
thetamin  =  0; 
thetamax  =  pi; 
deltatheta  =  (thetamax  -  thetamin)  /  (numtheta  -  1); 
theta  =  thetamin; deltatheta: thetamax;  %  angle  vector 

% - -FKXjBm  PREPAPATICN - 

rnaxn  =  15;  %  sum  to  maxn  rather  than  infinity 

ivec(0+l)=l; 

epsm(0+l)  =  1; 

minuslvec (0+1)  =  1; 

n2vec(0+l)  =  1; 

for  n=l:maxn 

ivec  (n+1)  =  i*ivec  (n) ; 
minuslvec  (n+1)  =  (“D'^n; 


%  nondimensional  wavenumber 
%  observation  points  are  at 
%  one  and  two  boat  lengths 


n2vec  (n+l)  =  2*n+l; 
epsm(n+l)  =  2; 
end  %  for  n — 1 

fact  (0+1)  =1; 
for  rf=l:2*maxn+l 

fact  (n+l)  =  n*fact(n); 
end  %  n — 2 

for  n=0:nBxn 
for  rrt=0:maxn 

leg_fact_rat  (n+l,m+l)  =0; 
coeff2  (n+l,in+l)  =0; 
coeff3  (n+l,mfl)  =0; 
coeff4  (n+l,m+l)=0; 
end  %  for  m — 1 
end  %  for  n — 3 
for  n=0:maxn 
for  nK):n 

leg_fact_rat  (n+l,m+l)  =fact  (n-m+l)  /fact  (n+m+1)  ; 
coeff2  (n+l,m+l)=2*  ( (n-m+l)  *fact  (n+m+1+1) )  /  ( (2*n+3)  *fact  (n-m+l) )  ; 
coeff3  (n+l,m+l)  =2*  ( (n+m)  *  (n-m)  *fact  (n+m+1) )  /  ( (2*n-l)  *fact  (n-m+l) ) 
coeff4  (n+l,m+l)=2*  ( (n-m)  *  (n-m-l)  *fact  (n+m+1) )  /  ( (2*n-l)  *fact  (n-m+l) ) 
end  %  for  m — 2 
end  %  for  n — 4 

coeff2  ( : , 0+1) =coeff2 ( : ,  0+1) /2; 
coeff3 ( : , 0+1) =coeff3 ( : , 0+1) /2; 
coeff 4 ( : , 0+1) =coef f 4 ( ; , 0+1) /2; 

for  n=0:raaxn+l 
for  j=0:itaxn+l 

POcoeff (n+l, j+1)  =0; 
end  %  for  j — 1 
end  %  for  n — 5 

for  n=0:inaxn+l 
for  j=n-l;-2:0 

POcoeff (n+l, j+l)=l/ ( j* ( j+l)-n* (n+l) ) ; 
end  %  for  j — 2 
end  %  for  n — 6 
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% - PPOGRAM - 


% - BESSEL  CALL - 

[jn  djn  yn  cfyn]  =  bessel_sph(maxn,kRK:a) ; 
hn  =  jn  +  i*yn; 
dhn  =  djn  +  i*cfyn; 

kappa_l  =  ivec. *  (djn  { : ,  3) .  *hn  ( : ,  1) . /dhn  ( : ,  3) ) . ' ; 
k^pa_2  =  ivec.*  (djn(:,3)  .*hn(:,2)  ./dhn(:,3) ) . 

TnTj_l  =  k^pa_l . '  *con  j  (kaFpa_l)  +con j  (kappa_l . ' )  *kaFpa_l; 

TnT  j_2  =  k^pa_2 . '  *con  j  (kaFpa_2)  +con  j  (kaFpa_2 . ' )  *kaFpa_2; 

TnTn_l  =  kaFpa_l.*conj  (kaFpa_l) ; 

TnTn_2  =  kappa_2.*conj  (kappa_2) ; 

[P0,dP0]  =  LegendrePnm(maxn+l,0) ; 

%  - SUMl - 

siird_l===n2vec*TnTn_l . '  ; 
smil_2==n2vec*TnTn_2 . ' ; 

%  — Theta  loop — 

for  k=l;numth€ta 

[P  dP]  =  LegendrePnni(maxn,  cos  (theta  (k) ) ) ; 

Pcx)s  =  leg_fact_rat.*P; 

%  - suyi2 - 

sum2_l=0; 
sum2_2=0; 
for  n=0:maxn 
STjn2iit=0; 
for  irH):n 

s\rn2rrF=sim2rn+Pcos  (n+l,m+l)  *Pcos  (n+l,iTH-l)  *coeff2  (n+l,m4-l) ; 
end  %for  m — siot2 
sura2_l===s\rri2_l+TnTn_l  (n+1)  *sum2m; 
sum2_2==Biim2_2+TnTn_2  (n+1)  *sum2m; 
end  %for  n — suin2 
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•SUM3- 


%  - 

sum3_l=0; 
sijm3_2=0; 
for  n=2:inaxn 
si]m3m=0; 
for  in=0:n“2 

siini3te=sim^Pcos  {n+l,rrH-l)  *Pcos  (n-2+l,m+l)  *coeff3  (n+l^m+l) ; 
end  %for  m — sumS 

sim6_l=siim3_l+TnTj_l  (n+l,n-2+l)  *si]m3m; 
surn3_2=surn3_2+TnTj_2  (n+l,n-2+l)  *suiii^; 
end  %for  n — sum3 

%  - SUM4 - 

sum4_l=0; 
sum4_2=0; 
for  n=l:inaxn 
sum4irf=0; 
for  iTi=0:n-l 

sum4m=sim4rrH-Pcos  (n+l,iTH-l)  *Pcos  (n+l,m+l)  *coeff4  (n+l,nM-l)  ; 
end  %for  m — sum4 
sum4_l=sum4_l+TnTn_l  (n+1)  *sijm4m; 
sum4_2==siin4__2+TnTn_2  (n+1)  *siim4m; 
end  %for  n — sum4 

% - SUMS - 

sum5_l=0; 
sum5_2=0; 
for  n  =  l:maxn 
for  j  =  n-l:”2:0 
sum5tn=0; 
for  iifO:  j 

POtem^PO  (n+l+l,in+l)  *dP0  { j+l+l,m+l)  -dPO  (n+l+l,iTH-l)  *P0  ( j+l+l,m+l)  ; 
ir!Coeff=epsm(m+l)  *  (n-m+l)  *  ( ; 

siin5ni=sum5rtH-Pcos  (n+l,m+l)  *Pcos  ( j+l,m+l)  *incoeff*POtenn; 
end  %for  m — s\jrr6 

sum5_l  =  sum5_l+si2m5m*TnTj_l  (n+1,  j+1)  *P0coeff  (n+1+1,  j+1+1)  ; 
sum5_2  =  sxjm5_2+siim5m*TnTj_2  (n+1,  j+1)  *P0coeff  (n+1+1,  j+1+1) ; 
end  %for  j — sumS 
end  %for  n — sumS 
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SUMS' 


%  - 

sum6_l=0; 
sum6_2=0; 
for  n  =  2:niaxn 
for  j  =  n-l:-2:l 
sum&KpO; 
for  m=0:  j“l 

P0term=P0  (n+l+l,m+l)  *dP0  ( j-l+l,m+l)-dP0  (n+l+l,m+l)  *P0  ( j-l+l,m+l)  ; 
mcoeff=^sm(irH-l)  *  (n-m+l)  *  ( j+m) ; 

sum6in=sum6m+Pcos  (n+l,m+l)  *Pcos  ( j+l,iTH-l)  *incoeff*POtem; 
end  %for  m — sumS 

sum6_l  =  sum6_l+sum6m*TnTj_l  (n+1,  j+1)  *P0coeff  (n+1+1,  j-1+1)  ; 
siim6_2  =  siin6_2+sum6m*TnTj_2  (n+l,  j+1)  *P0coeff  (n+1+1,  j-1+1) ; 
end  %for  j — sum6 
end  %for  n — sum6 

%  - SUM7 - 

sum7_l=0; 
sum7_2=0; 
for  n  =  Itmaxn 
for  j  =  n-l:-2:0 
sum7m=0; 
for  nt=0:  j 

P0term=P0  (n-l+l,m+l)  *dP0  ( j+l+l,m+l)  -dPO  (n-l+l,m+l)  *P0  ( j+l+l,iTH-l)  ; 
incoeff=epsni{irt+l)  *  (n-nn)  *  ( j-m+1) ; 

sum7m=siini7m+Pcos  {n+l,iT+'l)  *Pcos  ( j+l,nvfl)  *mcoeff*POterm; 
end  %for  m — sum7 

suni7_l  =  s'um7_l+sum7m*TnTj_l  (n+1,  j+1)  *P0coeff  (n-1+1,  j+1+1) ; 
surti7_2  ==  sum7_2+sum7m*TnTj_2  (n+1,  j+1)  *P0coeff  (n-1+1,  j+1+1) ; 
end  %for  j — sum7 
end  %for  n — sum7 

%  - SUMS - 

suin8_l=0; 
s\jn8_2=C; 
for  n  =  2:iiiaxn 
for  j  =  n-l:-2:l 
sumSnK); 
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for  m=0:  j-1 

P0tenTF=P0  {n-l+l,m+l)  *dP0  ( j-l+l,m+l)  -dPO  (n-l+l,m+l)  *P0  ( j-l+l,m+l)  ; 
mcoeff=epsm(rrH-l)  *  (n+m)  *  { j+m) ; 

sijmarrt=sum8m+Pcos  (n+l,n:H-l)  *Pcos  { j+l,m+l)  *mco€ff*POterm; 
end  %for  m — simS 

sum8_l  =  sum8_l+sum8m*TnTj_l  (n+1,  j+1)  *P0coeff  (n-1+1,  j-1+1) ; 
sum8_2  =  sum8_2+sum8m*TnTj_2(n+l,  j+l)*P0coeff  (p-1+1,  j-1+1); 
end  %for  j — si3m8 
end  %for  n — sum8 

pl_l=suml_l; 

pl_2=suirLL_2; 

p2_l=(surn2_l+suiii3_l+sijrn4_l)  ; 
p2_2=(sum2_2+sum3_2+sum4_2) ; 
p3_l=(surr5_l+sum6_l+si]m7_l+sum8_l) ; 
p3_2=(siir5_2+sum6_2+sum7_2+suin8_2) ; 

Pressurel  (k)  =4*pi*alpha*pl_l+2*pi*  (beta+gairma)  *p2_l+2*pi*  (beta-gainm)  *p3_l; 
Pressure2  (k)=4*pi*alpha*pl_2+2*pi*  (beta+ganma)  *p2_2+2*pi*  (beta-gainna)  *p3_2; 
end  %  for  theta 

pn=2*pi*  (2*alpha+  (1/3)  *  (beta+gamna) ) ; 

TSl=10*logl0  (Pressure/pn)  ; 

TS2=10*logl0  (Pressure2/pn) ; 

TS(1,:)=TS1; 

TS(2,:)=TS2; 

end 
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